Breakdown of the large-scale wind in F = 1/2 rotating Rayleigh-Benard flow 
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Experiments and simulations of rotating Rayleigh-Benard convection in cylindrical samples have 
revealed an increase in heat transport with increasing rotation rate. This heat transport enhance- 
ment is intimately related to a transition in the turbulent flow structure from a regime dominated 
by a large-scale circulation (LSC), consisting of a single convection roll, at no or weak rotation to 
a regime dominated by vertically-aligned vortices at strong rotation. For a sample with an aspect 
ratio r = D/L = l(Dis the sample diameter and L its height) the transition between the two 
regimes is indicated by a strong decrease in the LSC strength. In contrast, for F = 1/2 Weiss and 
Ahlers [J. Fluid Mech. 688, 461 (2011)] revealed the presence of a LSC-like sidewall temperature 
signature beyond the critical rotation rate. They suggested that this might be due to the formation 
of a two- vortex state, in which one vortex extends vertically from the bottom into the sample interior 
and brings up warm fluid, while another vortex brings down cold fluid from the top; this flow field 
would yield a sidewall temperature signature similar to that of the LSC. Here we show by direct 
numerical simulations for F = 1/2 and parameters that allow direct comparison with experiment 
that the spatial organization of the vertically-aligned vortical structures in the convection cell do 
indeed yield (for the time average) a sinusoidal variation of the temperature near the sidewall, as 
found in the experiment. This is also the essential and non-trivial difference with the F = 1 sample, 
where the vertically-aligned vortices are distributed randomly. 



I. INTRODUCTION 

The flow of a fluid heated from below and cooled from 
above, better known as Rayleigh-Benard (RB) convec- 
tion [H ^2 , is the classical system to study heat trans- 
fer phenomena. For given aspect ratio F = D/L (D is 
the sample diameter and L its height) and given geom- 
etry, its dynamics is determined by the Rayleigh num- 
ber Ra=/3^AI/^/(/^i^) and the Prandtl number Vi=v / k.. 
Here, /3 is the thermal expansion coefficient, g is the 
gravitational acceleration, A is the temperature differ- 
ence between the horizontal plates, and v and n are the 
kinematic viscosity and thermal diffusivity, respectively. 
The response of the system is expressed by the dimen- 
sionless heat transfer, that is the Nusselt number Nu, and 
the Reynolds number Re. We will first give a brief in- 
troduction into the concept of the large scale circulation 
(LSC) (Section [lA| ) before we provide a brief summary of 
some important concepts and properties of rotating RB 



convection in Section ffB) Subsequently, in Section IC 



discuss the interesting aspects that are addressed in the 
present paper. 



plumes with cold fluid flow down on the opposite side of 
the convection cell setting up a large-scale mean flow in 
the cell. This fluid motion is better known as the LSC, 
see the diagram in figure [l^. A recent review on the LSC, 
its properties and dynamics is provided in Ref. [1]. 
Keeping the sketch of the LSC in mind we may expect 
hot rising fluid on one side and cold sinking fluid at the 
opposite side of the cell which may be visualized by hor- 
izontal temperature snapshots. A visualization of an in- 
stantaneous temperature field at mid height obtained in 
a simulation indeed gives an impression of the structure 
of the LSC, see figure [TJd. In experiments the LSC is 
measured by thermistors in the sidewall, see for exam- 
ple Brown et al. [3 , or by small thermistors that are 
placed in the flow at different azimuthal positions and 
different heights, see for example Xia et al. [4 . Since the 
LSC transports warm (cold) fluid from the bottom (top) 
plate up (down) the sidewall, the thermistors can detect 
the location of the up-flow (down-flow) by showing a rel- 
atively high (low) temperature effectively resulting in a 
sinusoidal azimuthal temperature profile. 



A. The LSC in Rayleigh-Benard convection 

Due to the temperature difference between the hor- 
izontal plates warm plumes rise from the hot bottom 
plate and cold plumes sink from the cold top plate. In 
our cylindrical convection cell the collection of rising 
and sinking plumes organize as follows: the plumes with 
warm fluid collect and flow up on one side of the cell and 



In order to analyze the sidewall temperature profile 
Stevens et al. [6 proposed a quantitative measure for 
the temperature signature, which they called the relative 
LSC strength in their study on the effect of plumes on 
measuring the LSC in samples with F = 1 and F = 1/2 
in non-rotating RB convection. This measure is based on 
the energy in the different Fourier modes of the measured 
or computed azimuthal temperature profile at or nearby 



(a) (b) (c) 

FIG. 1. (Color online) a) Sketch of a single roll state (left) and a double roll state (right) of the large scale circulation in a F = 1/2 
sample. (b,c) Visualization of the temperature field in the horizontal mid-plane of a cylindrical convection cell for Ra=2.91 
xlO^ and Pr=4.38 in a F = 1 sample 5 . The red and blue areas indicate warm and cold fluid, respectively, (b) At a slow 
rotation rate (1/Ro=0.278), and below the transition: The warm up-flow (red) and cold down-flow (blue) reveal the existence 
of a single convection roll superimposed upon turbulent fluctuations, (c) At a somewhat larger rotation rate (1/Ro=0.775), 
above the transition: Here, the vertically- aligned vortices cause disorder on a smaller length scale in the horizontal plane. 



the sidewall, as 

5fc = maxM^^^,0 . (1) 

The subscript k indicates the height at which Sk is de- 
termined, that is for /c = 6 at z/L = 0.25, for /c = m at 
z/L = 0.50, and for k = t at z/L = 0.75. In equation ([T]) 
(El) = J^^ Eidt indicates the sum of the energy in the 
first Fourier mode over the time interval [t^, t^] (from the 
beginning of the statistically steady part of the simula- 
tion at t = ^5 to the end of the simulation at t = te)^ and 
{Efot) is defined similarly with Etot being the total energy 
in all Fourier modes. Furthermore, N is the total number 
of Fourier modes that can be determined (depending on 
the number of azimuthal probes). From its definition it 
follows that < Sk < I- Concerning the limiting values: 
Sk = I means the presence of a pure azimuthal cosine 
profile and Sk = indicates that the magnitude of the 
cosine mode is equal to (or weaker than) the value ex- 
pected from a random noise signal, see Ref. [6^. Hence, 
a value for Sk of about 0.5 or higher is a signature that 
a cosine fit on average is a reasonable approximation of 
the data set, as then most energy in the signal resides 
in the first Fourier mode. In contrast, Sk well below 0.5 
indicates that most energy resides in the higher Fourier 
modes. In section [TbI we discuss how this method is used 
to determine the existence of different turbulent states in 
rotating RB convection. 



B. Rotating Rayleigh-Benard convection 

The case where the RB system is rotated around its 
vertical axis at an angular speed ft is used to study the 
effect of rotation on heat transport and flow structuring. 
Here we define the Rossby number as Ro = U/{2QL), 
with U = ^/Pg/S.L being the free-fall velocity (mean bulk 
velocity) to make sure that Ro is the relevant parameter 
that determines when the formation of large scale vortices 
sets in. In the remainder of the paper we will indicate 
the dimensionless rotation rate by 1 /Ro as this parameter 
increases with increasing rotation rate. 

Since the experiments by Rossby [7 it is known that 
rotation can enhance heat transport. Rossby measured 
an increase in the heat transport with respect to the non- 
rotating case of about 10% when water is used as the 
convective fluid. This increase is counterintuitive as the 
stability analysis of Chandrasekhar [8 has shown that ro- 
tation delays the onset to convection and from this anal- 
ysis one would expect that the heat transport decreases. 
The mechanism responsible for this heat transport en- 
hancement is Ekman pumping |7( l9HT3j. that is, due to 
rotation, rising or falling plumes of hot or cold fluid are 
stretched into vertically-aligned vortices that suck fluid 
out of the thermal boundary layers adjacent to the bot- 
tom and top horizontal plates. Evidence for the existence 
of vertically- aligned vortices was reported 15-20 years ago 
p!4HT6] . Sakai was the first who confirmed with flow vi- 
sualization experiments that there is a typical ordering 
of vertically-aligned vortices under the influence of rota- 
tion. For higher rotation rates a strong heat transport 
reduction, due to the suppression of the vertical velocity 
fluctuations by rotation, is found. After Rossby [7] many 
experiments have confirmed this general picture quali- 
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tatively and, in recent years, also quantitatively, see for 
example Refs. ^ US UZHH] . 

Next to experiments there have been a number of numer- 
ical studies of rotating RB convection [F, ^ [HI HH |24]- 
[29j . These studies focussed on the influence of rotation 
on heat transport and the corresponding changes in the 
flow structure. Zhong et al. [13 and Stevens et at. j5l[27] 
used results from experiments and simulations in a F = 1 
sample to study the influence of Ra and Pr on the effect 
of Ekman pumping, and thus on heat transport. 
For an overview of the parameter regimes that are cov- 
ered in simulations and experiments during the recent 
decades we refer to the flgures 1 of Stevens et al. [28l [29] . 
These phase diagrams reveal that there are two ap- 
proaches to analyze rotating RB convection. The flrst 
approach is to vary the rotation rate while the tempera- 
ture difference between the two horizontal plates is flxed. 
This means that one varies 1/Ro, while Ra is kept con- 
stant. The second approach is to keep the ratio between 
the viscous force to the Coriolis force, indicated by the 
Ekman (Ek) or Taylor (Ta) number, flxed. In almost 
all simulations and experiments the Pr number is kept 
constant and a flxed Ek number then means that Ra and 
1/Ro are varied. This approach is followed by several au- 
thors [T2 l [T8 lfT9t [3QH33]. Here we take the flrst approach, 
that is, vary the rotation rate while the temperature dif- 
ference between the plates is flxed. 

When heat transport enhancement as function of the 
rotation rate is considered, a typical division in three 
regimes is observed. Namely regime I (weak rotation), 
where no heat transport enhancement is observed, regime 
II (moderate rotation), where a strong heat transport 
enhancement is found, and regime III (strong rotation), 
where the heat transport starts to decrease [HI [24l [25] . 
Based on an experimental and numerical study on the 
properties of the LSC in rotating RB in a F = 1 sample 
Kunnen et al. [TT] showed that there is no LSC in regimes 
II and III. Subsequently, Stevens et al. [5 showed for a 
similar setting that the heat transport enhancement at 
the start of regime II sets in as a sharp transition (at 
a critical value of the inverse Rossby number, 1/Roc). 
They showed that in experiments the transition is indi- 
cated by changes in the time-averaged LSC amplitudes, 
i.e. the average amplitude of the cosine flt to the az- 
imuthal temperature proflle at the sidewall, and the ver- 
tical temperature gradient at the sidewall. We note that 
the formation of a mean vertical temperature gradient 
in rotating RB convection was already observed earlier, 
see Refs. [171 Later experiments and simulations 

[23[ [24] in a F = 1 sample revealed that Sk is close to 
one before the onset of heat transport enhancement. In 
this regime the LSC is the dominant flow structure. Af- 
ter the onset vertically-aligned vortices become the domi- 
nant flow structure and then Sk quickly decreases to zero. 
The reason is the formation of many randomly positioned 
vertically- aligned vortices, and due to their random loca- 
tions, see flgure [ij^, the cosine mode in the azimuthal 
temperature proflle disappears. 



In addition, simulation results have shown that the 
transition between the two different states is visible not 
only in the Nu number and the LSC statistics, but also 
by a strong increase in the number of vortices at the 
thermal boundary layer height [S] [2TJ [26l [J9 . In addi- 
tion, rotation changes the character of the kinetic bound- 
ary layer from a Prandtl-Blasius boundary layer at no or 
weak rotation to an Ekman boundary layer after the on- 
set, which is revealed by the 1/Ro ' scaling of the ki- 
netic boundary layer thickness after the onset for F = 1. 
Furthermore, Stevens et al. [5 observed for this case an 
increase in the vertical velocity fluctuations at the edge of 
the thermal boundary layer, due to the effect of Ekman 
pumping, while a decrease in the volume averaged verti- 
cal velocity fluctuations is found due to the destruction 
of the LSC. 



C. Issues addressed in the present paper 

Most flndings described above are based on experi- 
ments and simulations in F = 1 samples. In the present 
paper we study rotating RB convection in an aspect ratio 
F = 1/2 sample as experiments have revealed important 
differences with respect to the F = 1 case. Recently, 
Weiss and Ahlers [35 revealed that in such a F = 1/2 
sample with water (Pr=4.38) at a range of Rayleigh num- 
bers (2.3 X 10^ <Ra< 7.2 x 10^°) the flow can be either 
in a single roll (see left panel of flgure [l^) or in a double 
roll state (see right panel of flgure [l^) when no rotation 
is applied. Here, modest rotation (l/Ro< 1/Roc), stabi- 
lizes the single roll state and suppresses the double roll 
state. Surprisingly, computation of Sk reveals that after 
the onset of heat transport enhancement at 1/Roc the 
cosine-signature of the temperature proflle at the side- 
wall does not disappear. This may suggest that in a 
F = 1/2 sample the single roll state continues to exist in 
the rotating regime, whereas it has been shown to disap- 
pear in a F = 1 sample [S] [TTJ |23l [24] . As discussed by 
Weiss and Ahlers [H [23] the value of 5/^ > 0.5 does not 
allow us to distinguish between a single roll state and a 
two- vortex state, in which one vortex extends vertically 
from the bottom into the sample interior and brings up 
warm fluid, while another vortex transports cold fluid 
downwards at the opposite part of the cylindrical sam- 
ple. Such a two- vortex state presumably results in a pe- 
riodic azimuthal temperature variation close to the side- 
wall which cannot be distinguished from the temperature 
signature of a convection roll with up-flow and down-flow 
near the side wall. It is worthwhile to emphasize the dis- 
tinction with the F = 1 case with l/Ro>l/Roc. Here, the 
vertically-aligned vortices emerging at the bottom plate 
(transporting hot fluid upwards) are distributed more or 
less randomly over the horizontal cross section and simi- 
larly for the vortex tubes emanating from the top plate. 
From a retrospective point of view we can now conclude 
that Stevens et al. [29] already observed some flrst hints 
for the presence of a two- vortex state (see flgure [s] for an 
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TABLE I. Details for the simulations performed at Ra=4.52 x 10^ and Pr=4.38 in an aspect ratio F = 1/2 sample. The columns 
from left to right indicate the Rossby number 1/Ro, the Nusselt number (Nuf) obtained after averaging the results of the 
three methods (see text) using the whole simulation length, and the Nusselt number (Nuh) after averaging the results of the 
three methods over the last half of the simulation. The following column indicates the normalized temperature gradient at the 
sidewall Aw/ A. The next three columns give the Nusselt number derived from the volume averaged value of the pseudo kinetic 
dissipation rate (e^^), and the kinetic (cu) and thermal (ee) dissipation rates compared to Nuf indicated in column two. The 
last column indicates the simulation length in dimensionless time units that is considered for the data collection; before this 
period each simulation has been run for 80 — 100 time units in order to prevent transient effects. 
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impression of a two- vortex state) in a F = 1/2 sample at 
Ra = 2.91 X 10^ 1/Ro = 3.33 and Pr = 4.38 (see figure 
5a of that paper). However, an open question is whether 
this state prevails for a considerable time and for differ- 
ent values of Ra, Pr, and a range of rotation rates (in 
particular matching the settings of the experiments by 
Weiss and Ahlers [22l [23| ) . Additionally, can we comple- 
ment the indirect evidence for a two- vortex state as put 
forward by Weiss and Ahlers [22] |23] with direct obser- 
vations of the flow structures in the bulk flow? Can we 
also distinguish between a genuine two- vortex state and 
a flow structuring consisting of a few vertically- aligned 
vortices with hot rising fluid on one side of the cylindri- 
cal sample and a few of such vortices bringing cold fluid 
from the top plate downwards close to the opposite side- 
wall? In the present paper we report on direct numerical 
simulation results for F = 1/2, Ra = 4.52 x 10^, and Pr 
= 4.38 over the range < 1/Ro < 12, corresponding to 
a subset of the Weiss-Ahlers experiments thus allowing 
direct comparison. The results show that a complex vor- 
tex state, which in the time average has the signature of 
a two- vortex state, persists over the range 1/Roc ^ 1/Ro 
< 12. This confirms that the observation of a sinusoidal 
azimuthal temperature profile near the sidewall can in- 
deed be explained by the presence of Ekman vortices. In 
addition, to this general result, we report detailed data 
that permit a direct comparison of 5^, Sm^ and 5*5 (see 
figure [3| and the LSC amplitude (see figure [4| with ex- 



periment over a wide 1/Ro range. Furthermore, we also 
computed the Nusselt number and the vertical tempera- 
ture gradient at the sidewall (see figure |2|. All of these 
properties show excellent agreement between our simula- 
tions and the available experimental measurements; our 
simulations can therefore be used for exploration of the 
bulk and boundary layer flow structure which are not 
accessible in the current experiments. 

We first discuss the numerical method in Section |n| 
before we compare the simulation data with the exper- 
iments in Section IIIII In Section IIVI we discuss flow di- 
agnostics that can only be obtained in simulations and 
show that in a F = 1/2 sample the vertically-aligned 
vortices arrange such that a sinusoidal azimuthal tem- 
perature profile close to the sidewall is formed. 



II. NUMERICAL METHOD 

We simulate rotating RB convection for Ra=4.52 x 10^ 
and Pr=4.38 in an aspect ratio F = 1/2 sample by solv- 
ing the three-dimensional Navier- Stokes equations within 
the Boussinesq approximation. A constant temperature 
boundary condition is applied at the horizontal plates, 
while the sidewall is modeled as adiabatic. This case has 
been chosen to be as close as possible to the experiments 
performed by Weiss and Ahlers [22j |23]. General de- 
tails about the numerical procedure can be found in Refs. 
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[36l [37] and specific details concerning the (non) rotating 
RB simulations in Refs. [ 13l |27]. 

In order to eliminate the effect of transients we dis- 
carded the information of the first 80 — 100 dimension- 
less time units and the simulation lengths we mention 
refer to the length of the actual simulation, thus the 
period after this initialization. In Table [T| we compare 
the Nusselt number averaged over the whole simulation 
length (denoted by Nuf] the Nusselt number is based on 
the averages of three methods, i.e. the volume average 
of Nu = {{uzO)a — i^AL~^ and the averages 

based on the temperature gradients at the bottom and 
top plate) with the Nusselt number averaged over half 
the simulation length {Nuh)- For all cases these values 
Nuf and Nuh are converged within 1%. 
The simulations have been performed on a grid with 
641 X 161 X 641 nodes in the azimuthal, radial, and axial 
direction, respectively. The grid allows for a very good 
resolution of the small scales both inside the bulk of tur- 
bulence and in the boundary layers where the grid-point 
density has been enhanced. We checked this by calculat- 
ing the Nusselt from the volume- averaged kinetic energy 
dissipation rate (e^) = v^{Nu — l)Ra/{Pr'^L'^)^ and ther- 
mal dissipation rate (e^i) = n/S?Nu/ L'^ as is proposed by 
Stevens et al. [38]. In addition, we now also compare 
with the volume averaged value of := u • V^u. The 
Nusselt number calculated from these quantities is al- 
ways within a 5% margin, and even much closer for most 
simulations, of Nuf^ which according to Stevens et al. 
[38 indicates that the simulations are well-resolved, see 
Table [l| for details. As argued by Shishkina et al. [39 it 
is especially important to properly resolve the boundary 
layers. Our grid-point resolution in the boundary lay- 
ers also satisfy their criteria for the rotating case (where 
kinetic boundary layers tend to become thinner with in- 
creasing rotation rate). 



III. COMPARISON WITH EXPERIMENTS 

In figure [2^ we show that the Nusselt number as func- 
tion of the rotation rate Nu(l/Ro), with respect to the 
non-rotating value Nu(0), obtained in the simulations 
agrees excellently with the experimental result of Weiss 
and Ahlers [23]. In addition, we note that the absolute 
values differ less than 1%, which can be considered as 
an excellent agreement. The figure shows that a strong 
heat transport enhancement due to Ekman pumping sets 
in at a critical dimensionless rotation rate 1/Roc ~ 0.86. 
For strong rotation rates, i.e. high values of 1/Ro, the 
expected decrease in the heat transport is observed. In 
figure |2]3, we find that the normalized vertical temper- 
ature gradient at the sidewall at z/L = 0.50, denoted 
by A^y/A, and calculated from the azimuthally and time 
averaged temperatures at the sidewall at z/L = 0.25 and 
z/L = 0.75, is also in good agreement with the measure- 
ments by Weiss and Ahlers [23 . 

In the simulations we placed 64 numerical probes at 
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FIG. 2. (color online) (a) The heat transfer enhancement 
with respect to the non-rotating case, Nu(l/Ro)/Nu(0), and 
(b) the temperature gradient at the sidewall as function of the 
rotation rate 1/Ro for Ra=4.52 x 10^ and Pr=4.38 in a T = 
1/2 sample. The experimental data from Weiss and Ahlers 
[23] are indicated by the blue open squares and the simulation 
results with the red solid circles. The vertical dashed lines 
indicates the position of the transition (at 1/Roc ~ 0.86) [2T] . 



r = 0.951^/2, which is outside the sidewall boundary 
layer. The probes are placed some distance from the 
wall in order to collect data on the velocity components 
as well, which are all zero at the wall. These velocity data 
are complementary to the temperature measurements as 
the region in which a relatively high (low) temperature is 
measured should correspond to the region with positive 
(negative) vertical velocity. As we do not want to see 
the effect of very small plume events we apply a moving 
averaging filter of 8 dimensionless time units, see Stevens 
et al. [6] for details, to the temperature measurements of 
the probes, before we determined Sk, with k G {6, m,t}. 
Figure [3^-c shows that the measured Sk at z/L = 0.25, 
z/L = 0.50, and z/L = 0.75 in the simulations agree 
well with the experimental measurements [23J , which are 
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FIG. 3. (color online) The values of (a) Sb (at 0.25z/L), (b) Sm (at 0.50z/L), and (c) St (at 0.75z/L) as function of 1/Ro. 
The experimental data from Weiss and Ahlers [23 are indicated by blue open squares and the simulation results with red solid 
circles. The experimental and simulation results are for Ra=4.52 x 10^ and Pr=4.38 in a F = 1/2 sample. For the simulation 
data S is based on the temperature measurements of 64 azimuthally equally spaced probes outside the sidewall boundary layer 
and for the experimental data it is based on temperature measurements of 8 probes embedded in the sidewall. The vertical 
dashed line indicates the position of 1/Roc ~ 0.86 [21_. 



0.035 




1/Ro 



FIG. 4. (color online) The time- averaged temperature ampli- 
tudes (Sk) for all three levels k G {6, m,t}, with the symbols 
representing b (red circles) , m (black squares) and t (blue tri- 
angles)} as function of 1/Ro for F = 1/2 and Ra=4.52 x 10^. 
Vertical dashed line as in figure |3] 



numerical probes at the height k G {6, m,t}, with 6, m, 
and t defined below Eq. ([T]). Here, (pi = 2z7r/A/', with N 
the number of probes, refers to the azimuthal position 
of the probes, and and (j)^ indicate the temperature 
amplitude and orientation of the LSC, respectively. In 
figure [4] we plot the time-averaged temperature ampli- 
tude {Sk) of the LSC as function of the rotation rate 
for r = 1/2. A comparison with the experimental data 
of Weiss and Ahlers [23] (where results are available for 
both F = 1/2 and F = 1) shows that within our statisti- 
cal convergence the trends shown in the numerical data 
for F = 1/2 are similar to those revealed in the experi- 
ments (and show a remarkably different trend compared 
to results for F = 1 which show a decreasing (Sk) for 
increasing rotation [23 ). The temperature amplitude of 
the LSC is relatively small for weak rotation (small val- 
ues of 1/Ro). It becomes slightly larger when increasing 
the rotation rate from zero to the critical rotation rate 
(1/Roc ^ 0.86). Subsequently, a small dip just after the 
onset of heat transport enhancement is observed, which 
is followed by a small further increase in the temperature 
amplitude. 



based on the temperature measurements of just 8 probes 
embedded in the sidewall instead of the 64 numerical 
probes in our simulations. The error bars in the figure in- 
dicate the difference between Sk obtained using the com- 
plete time interval with Sk based on the last half of the 
simulation. Experimental and numerical data obtained 
in a F = 1 sample revealed that Sk strongly decreases 
to values around Sk ^ 0.2 when the heat transport en- 
hancement sets in ^23j. This is obviously not the case in 
the F = 1/2 sample (in general Sk ^ 0.5 for 1 < 1/Ro < 
10, see figure[3|. The experimental results do not suggest 
that this general feature depends on Ra [23] . 

Following Ref. [40 the orientation and strength of the 
LSC can be determined by fitting the function Oi = 
Ok + Sk cos{(j)i — (pk) to the temperatures recorded by the 



In summary, in a F = 1 sample the transition from 
the regime with no or weak rotation, with the LSC as 
the dominant flow structure, to the rotation dominated 
regime, with vertically-aligned vortices, is indicated by 
a strong reduction of Sk- Also {Sk) decreases with in- 
creasing rotation rate also suggesting disappearance of 
the LSC. In contrast, in a F = 1/2 sample these criteria 
do not provide evidence that the LSC is destroyed at the 
onset of heat transport enhancement. In the following 
section we use data that are only available in the simu- 
lations to show that also in a F = 1/2 sample the LSC 
is destroyed at the onset of heat transport enhancement 
and that the vortices arrange such that on average a si- 
nusoidal temperature profile close to the sidewall is still 
measured. 




FIG. 5. (color online) Flow visualization for Ra=4.52 x 10^, 
Pr=4.38, and l/Ro=3.33 in an aspect ratio V — 1/2 sam- 
ple. The colormap for the horizontal planes at z/L — 0.25, 
z/L — 0.50, and z/L — 0.75 is different in panel a and b, 
and is indicated on the left and right hand side, respectively. 
The red temperature isosurface, originating from the bottom, 
indicates the region A > 0.65, whereas the blue isosurface, 
originating from the top, indicates the region A < 0.35. Cor- 
responding movies can be found in the supplementary mate- 
rial 1421. 



IV. FLOW STRUCTURES AFTER 

TRANSITION TOWARDS 
ROTATION-DOMINATED REGIME 



In this Section we use the availability of all flow data 
from the simulations to explain why in the rotating 
regime in an aspect ratio F = 1/2 sample the value of 
Sk does not decrease to values as small as 0.2 or less as 
usually observed in the F = 1 case. In particular, we flnd 
that the minimum value is Sk ^ 0.5 for 1 < 1/Ro < 10 
(see flgure[3|. Only for 1/Ro > 10 the experimental data 
indicate a further substantial decrease of Sk [23 . Weiss 
and Ahlers already proposed the idea of a two-vortex 
state which may explain the temperature proflles mea- 
sured in their experiments. However, no experimental 
validation was possible in their experiments. To address 
this issue we have visualized the flow and temperature 
fleld for Ra = 4.52 x 10^ and Pr=4.38 at 1/Ro = 3.33, 
see flgures[5](isosurfaces with constant temperature) and 
[g] (visualization of vortices by means of temperature and 
vorticity) and the complementary movies. The visualiza- 
tion based on temperature isosurfaces reveals that also 
at this high Ra number vertically- aligned vortices are 
formed in the rotation dominated regime. In the cor- 
responding movies one can see that close to the bottom 
(top) plate several vortices containing warm (cold) fluid 
are formed. Although several vortices are created close 
to the horizontal plates only a few vortices are strong 







FIG. 6. (color online) Flow visualization for Ra=4.52 x 10^, 
Pr=4.38, and l/Ro=3.33 in an aspect ratio F = 1/2 sample 
The top row shows three-dimensional temperature isosurfaces 
for two different time instances where the red region indicates 
A > 0.65, whereas the blue isosurfaces indicate the region 
A < 0.35, see figure [s] The bottom row shows a visualization 
of the vortices, based on the Q2,d criterion. Note that the red 
(blue) regions indicated by the temperature isosurface indeed 
correspond to vortex regions. 



enough to reach the horizontal planes (at z/L = 0.25, 
z/L = 0.50 and z/L = 0.75) so that they can be ob- 
served by sidewall temperature probes. Therefore, most 
of the time we see only a few prevailing vortices at the 
horizontal plane at midheight. In addition, we notice the 
non-trivial result that the vortices arrange such that all 
warm rising vortices tend to cluster on one side of the 
cell and the sinking cold vortices on the opposite side 
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(contrary to what is observed for F = 1). Note that 
the horizontal planes (at z/L = 0.25, z/L = 0.50 and 
z/L = 0.75) reveal that the main rising and sinking vor- 
tices are surrounded by a warm (cold) region. As these 
warm and cold regions are formed on opposite sides of 
the cell this explains the cosine temperature profile that 
is observed close to the sidewall. 

To further confirm that indeed vertically- aligned vortices 
are formed we analyzed the flow structure in more de- 
tail by determining the positions of the vortices by ap- 
plying a vortex detection algorithm. The fully three- 
dimensional detection technique is based on the velocity 
gradient tensor Vu = diUj{i^j G 1,2,3). This tensor 
can be split into a symmetric and antisymmetric part 
Vu = ^[Vu + (Vu)^] + ^[Vu - (Vu)^] = S^n. The 
QsD criterion, according to Hunt et al. [41 , defines a re- 
gion as vortex when Q2,d = ^(II^IP — II^^IP) > 0? where 
1 1^4 1 1 = ^JTr{AA^) represents the Euclidean norm of the 
tensor A. In practice the threshold value of zero to dis- 
tinguish between vorticity and strain dominated regions 
in the flow, and thus identifying vortices, is found to be 
unsuitable [26 . Here, we apply a relatively high posi- 
tive threshold value, identical for all presented snapshots, 
that is 'matched' to the temperature isosurfaces that are 
presented in figure [g] (top row) . The position and size of 
the detected vortex tubes also agree very well with the 
coherent structures indicated by the temperature isosur- 
faces, see figure |6] This confirms that the warm and cold 
regions shown by temperature isosurfaces are indeed vor- 
tices. The main difference between the vortices identified 
with the vortex detection algorithms and the areas indi- 
cated by the temperature isosurfaces are the smaller vor- 
tices in the middle. These small vortices do not show up 
when the temperature criterion is used, since their base 
is not close to the bottom (top) plate where warm (cold) 
fluid enters the vortices at their base. 
We computed the boundary-layer thickness as function of 
rotation rate and determined averaged root-mean-square 
velocity fluctuations, as has been done previously for the 
r = 1 case [t51 128] . Our aim is to confirm directly from 
flow field data that the flow structure in a F = 1/2 sam- 
ple indeed changes from a regime dominated by the LSC 
(1/Ro < 1/Roc) to one dominated by vertically-aligned 
vortices (1/Ro > 1/Roc). 

First of all we determined the kinetic boundary-layer 
thickness. Like in Ref. [28^ for the F = 1 sample the ki- 
netic boundary-layer thickness is approximately constant 
before the onset of heat transport enhancement sets in 
at the critical rotation rate 1/Roc. After the onset the 
kinetic boundary- layer thickness scales as (l/Ro)~^"^^, 
which is in agreement with the scaling expected from 
Ekman boundary layer theory. The current data for the 
F = 1/2 (not shown here; graph similar as for F = 1 [5j) 
thus reveal that the boundary-layer structure changes 
from a Prandtl-Blasius type boundary layer when no or 
weak rotation is applied to an Ekman type boundary 
layer when l/Ro>l/Roc. Subsequently, we determined 
the volume- averaged vertical velocity fluctuations, and 



the vertical velocity fluctuations at the edge of the ther- 
mal boundary layer. Although not shown here, our com- 
putations for F = 1/2 indicate that the volume-averaged 
vertical velocity fluctuations slightly increase for 1/Ro 
< 1/Roc and strongly decrease for 1/Ro > 1/Roc, thus 
indicating that the LSC is destroyed. 

The strong decrease in the normalized volume- 
averaged vertical velocity fluctuations coincides with a 
further increase of the horizontal average at the edge of 
the thermal boundary layers. The increase of the fluctua- 
tions at the boundary-layer height signifies enhanced Ek- 
man transport. Thus these averages provide additional 
support that the dominant flow structure changes from a 
LSC to a regime dominated by vertically-aligned vortices. 
The average arrangement of these few vortices supports 
the presence of a sinusoidal azimuthal temperature profile 
at (or close to) the sidewall for F = 1/2. This is also the 
essential and non-trivial difference with the F = 1 sam- 
ple, where the vertically- aligned vortices are distributed 
randomly. 



V. CONCLUSION 

We have shown that for F = 1/2 and 1/Ro > 0.86 there 
is no single-roll LSC. Instead, we found a complex state 
of vortices which, in the time average, yields a sinusoidal 
azimuthal temperature variation observed in the experi- 
ment over a wide range of rotation rates. We have com- 
pared data from direct numerical simulations performed 
at Ra=4.52 x 10^ with Pr=4.38 (water) in an aspect ratio 
F = 1/2 sample at different rotation rates with the ex- 
perimental results of Weiss and Ahlers ^23j . We find very 
close agreement in global properties, i.e. the measured 
Nusselt number, and local measurements, i.e. the verti- 
cal temperature gradient at the sidewall and the behavior 
of Sk{l/Ro). In contrast to the F = 1 case both Sk and 
the temperature amplitude of the LSC do not indicate 
any significant changes at the moment that heat trans- 
port enhancement sets in. Weiss and Ahlers [23 already 
discussed that Sk cannot distinguish between a single roll 
state and a two- vortex state, in which one vortex extends 
vertically from the top into the sample interior and brings 
down cold fluid, while another emanates from the bottom 
and introduces warm fluid. Here, we resolved this issue 
and show with a visualization of temperature isosurfaces 
and with the use of advanced vortex detection algorithms 
that at high Ra the formation of a flow structure with 
basically a few dominant vortices leads to a sinusoidal- 
like azimuthal temperature profile close to the sidewall. 
This state gives a high value of Sk at the three measure- 
ment heights due to the warm (cold) fluid that spreads 
out in the horizontal direction (due to horizontal trans- 
port of heat, see Stevens et al. [28]). This smoothens 
the temperature peak in the azimuthal temperature pro- 
file at the sidewall. The observation that Sk ^ 0.5 for 
very large rotation rates (up to l/Ro~ 10) in a F = 1/2 
sample indicates that the hot and cold vortices on aver- 
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age must align themselves such that the upgoing (warm) 
and downgoing (cold) vortices are on opposite sides of 
the cell, but of course a different organization can be 
formed at certain time instances. This organization of 
the vertically-aligned vortices in the F = 1/2 sample is 
a non-trivial difference with the F = 1 case. There, the 
vertically-aligned vortices are distributed randomly. We 
are not sure what physical mechanism causes this dif- 
ference between a F = 1 and F = 1/2 sample. At the 
moment we are considering stereoscopic particle image 
velocimetry measurements in Eindhoven in order to di- 
rectly visualize the two- vortex state in a F = 1/2 RB 
sample and to obtain statistics over a much longer time 
domain than in the simulations. We note that the visu- 
alization was not possible in the earlier experiments of 
Weiss and Ahlers [23 as these experiments focused on 
getting an accurate measurement of the heat transport 
and hope that these results can answer this question. 

In experiments Weiss and Ahlers [23 showed no sig- 
nificant difference between the (Ei) / (Etot) as function of 
1/Ro (see figure 5 of their paper) for different Ra. Based 
on this observation we do not expect a strong Ra number 
dependence of the S{l/Ro) curve at a given aspect ratio. 
At the moment there are no measurements or simula- 
tions available that studied the Pr number dependence 
of S{l/Ro) and new measurements of numerical simula- 
tions would be necessary to determine whether there is a 
strong Pr number effect. 



From previous investigations on rotating RB convec- 
tion in F = 1 cells it was concluded that Sk is a good 
indicator for the presence of the LSC. From the results 
reported here we conclude that this quantity cannot pro- 
vide a unique answer whether the LSC is present or not in 
turbulent rotating convection in F = 1/2 samples. The 
reason for this is that the vortices in a F = 1/2 align 
in such a way that on average a cosine like temperature 
profile is formed in the azimuthal direction along the side- 
wall. 
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